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ABSTRACT 


A design model for a once-through waste heat recovery 
unit with a segmented fin-tube arrangement was developed 
along with a simple model of a combined gas and steam 
(COGAS) turbine propulsion system. These models were inte- 
grated and applied in a computer program written in FORTRAN 
IV for the IBM 360-67 computer. Waste heat recovery unit 
designs were produced and tested at off-design conditions. 
Using the space constraints and power requirements of a 
Navy destroyer-type ship, one design was selected and 
employed to make estimates of possible fuel savings to be 


realized through the application of a COGAS system. 
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I. INTRODUCTION 


A. BACKGROUND 

The U. S. Navy has made a commitment to gas turbine 
propulsion systems, with two classes of gas turbine ships 
entering the fleet (DD-963 and FFG-7) and one class pro- 
posed (DDG-47) at the time of this writing. The DD-963 class 
has two shafts with two General Electric LM 2500 gas tur- 
bines powering each shaft. The FFG-7 class is a single- 
shaft ship with two LM 2500 gas turbines. The DDG-47 class 
moopulsion plant, as proposed, will..be.similar to. that of 
the DD-963 class. 

Currently, means for conserving fuel are being sought 
in all sectors of the economy and government. In most of 
our steam turbine-powered ships fuel conservation must be 
pursued largely through operating practices such as reduced 
steaming time or reduced speed when steaming. Several 
alternatives exist, however, for fuel consumption reduction 
for the gas turbine propulsion system. These alternatives 
involve combining another type of propulsion system with 
tne gas turbine in order to reduce the overall fuel consump- 
tion by either operating the combined systems in parallel 
or by operating the less expensive system in the cruise 
mode. Some examples of these combinations are: 

Gas Turbine and Diesel (CODOG, GODAG) 
Gas Turbine and Steam (COGAS) (steam system not 


separately fired) 
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Gas Turbine and Steam (COSAG) (steam system 
separately fired) 
The COmbined Gas turbine And Steam turbine (COGAS) system 
is probably the best of these alternatives with respect to 
fuel consumption, initial cost, maintenance, and manning 
requirements. 

The COGAS ship propulsion plant uses the exhaust heat 
from a gas turbine to generate steam in a waste heat re- 
covery unit (WHRU). This steam is used to drive a steam 
turbine which is operated in parallel with the gas turbine, 
for ship propulsion. A common reduction gear connects both 
turbines to the propeller shaft. Figure 1 illustrates the 
COGAS system. 

Several studies have examined the potential of the COGAS 
system for improving propulsion plant fuel efficiency. 
Giblon and Rolih [Ref. 1] studied the commercial ship 
application of COGAS. References 2 and 3 investigated the 
feasibility of the COGAS system for military ship application. 
In all three studies the COGAS system was designed around 
a gas turbine brake horsepower of about 23,000. Additionally, 
all three employed a drum-type boiler as the waste heat 
recovery unit. 

In this study, COGAS system designs were considered for 
a destroyer-sized ship with two shafts and two LM-2500 gas 
turbines for each shaft. One gas turbine was equipped with 
a WHRU. In the COGAS mode of operation the WHRU-equipped 


gas turbine was operated in parallel with the steam turbine, 
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powering one shaft. The other shaft was allowed to drag. 

It was assumed that the principal employment of the COGAS 
system would be for "cruise" conditions. That is, the COGAS 
system would be operated primarily for steady steaming. The 
horsepower range over which the COGAS system could be 
operated is, theoretically, limited only to the range of 
powers over which the gas turbine can be operated. Thus, 
the power range would vary from the COGAS system power 
output with the gas turbine at idle to the system power 
output attained with the gas turbine operating at its maxi- 
mum continuous power rating. 

When the COGAS system operating range is selected, the 
desired goal of fuel conservation must be balanced against 
system reliability and routine maintenance requirements. 

If only one gas turbine is equipped for COGAS operations, 
the system downtime will increase as the turbine is operated 
at higher power settings. In order to increase the time 
between major gas turbine maintenance requirements, two 
COGAS systems could be installed in a twin-shaft ship. The 
increased overall reliability which this redundancy would 
afford, however, would be at least partially offset by the 
higher initial cost, the increased manning requirements, 

and the added weight of two systems. 

A once-through counter-cross-flow heat exchanger was 
selected as the WHRU. Reference 6 states some of the 
principal advantages of the once-through steam generator 


over the drum-type boiler. These advantages are as follows: 
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lower initial cost, faster response, more compactness, 
and lower operating cost. Additiorfally, the once-through 
unit provides a somewhat simpler model than the drum-type 
heat exchanger. A schematic representation of a once- 
through WHRU installed in the exhaust ducting of a gas 


turbine is shown in figure 2. 


B. OBJECTIVES 

There were five main objectives for this thesis. 
1. Develop a design model for the once-through WHRU. This 
model was to size a WHRU and predict is performance given 
initial conditions such as feedwater inlet temperature, 
steam outlet temperature, WHRU operating pressure, inlet 
and outlet gas conditions, and gas flow rate. 
2. Consider several sets of initial conditions and design 
constraints (e.g., maximum frontal dimensions) and produce 
WHRU designs for each set. 
3. Evaluate these designs and integrate them with a simple 
COGAS system model in order to test the design model and 
to develop a clearer understanding of the WHRU performance 
in the COGAS system. 
4. From the evaluated designs, predict the performance 
of the COGAS system. 
5. Provide a framework for future studies involving a 
detailed, complete COGAS system cycle analysis and 


optimization. 


2. 
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II. MODEL DESCRIPTION 


A. OVERVIEW AND INITIAL CALCULATIONS 
The waste heat recovery unit (WHRU) is modeled as a 
once-through counter-cross-flow heat exchanger. The model 
is applied in a computer program for use in heat exchanger 
design for a combined gas and steam (COGAS) turbine pro- 
pulsion plant system. The same model could be adapted 
for use in a complete cycle analysis of the COGAS system. 
The WHRU model is divided into three principal sections: 
heating, boiling and superheating. The specified initial 
conditions for the model are: (1) gas temperature into 


the heat exchanger (Poi) (2) minimum average gas exit 


temperature from the heat exchanger (Tg4), (3) gas flow 
rate (m,), (4) water inlet temperature (Te5), (5) water/ 
steam pressure (P.)- Useful background material for the 


model formulation was obtained from references 1-5. 

Two simplifications should be noted at this point. First, 
the water-side pressure drop in the WHRU is neglected. 
Second, gas temperatures, evaluated at several points in 
the model, are average temperatures. Specifically, the 
average inlet and outlet gas temperatures for the heating 
and superheating sections and the average inlet and outlet 
gas temperatures for each pass in the boiling section are 
evaluated. In each of these calculations the gas tempera- 


ture distribution is neglected. Clearly, even if the gas 
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turbine exhaust flow is uniform in temperature as it 
approaches the WHRU, the gas temperature will not be uni- 
form after the first WHRU pass is encountered. Due to the 
varying gas-to-fluid temperature difference across each 
pass the gas temperature will be distributed non-uniformly 
after the first WHRU pass. One can, however, calculate the 
average gas temperature approaching a particular pass 

based on the heat release from the gas in the previous pass 
and the average gas temperature entering the previous pass. 
Averaged gas temperatures calculated in this manner are 
used throughout the model. 

Prior to the heating section, some initial calculations 
are performed to determine the water/steam mass flow rate 
mma. co Set the pinch point. In order to calculate the 
water/steam mass flow rate an energy balance is performed, 
using the starting conditions, as follows: 


( ) 


=ein TC T  -T 
2 S aps gl g4 


where C g is evaluated at the average gas temperature 


Tyavg = ete) gas? = The fluid mass flow rate 1s 


calculated from 


where Ney and hey are inlet and outlet fluid enthalpy 
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respectively. Intermediate fluid and gas temperatures are 
calculated as follows: 
1. The heating section is assumed to heat the water from 


inlet conditions to the saturation temperature, A 


Teo 
heat balance is performed on the water side yielding a 


heating section heat transfer of 


The required average gas temperature entering the heating 


section is then calculated from 


2. In the boiling section it is assumed that sufficient 
heat is added to produce saturated vapor at the outlet. 
An energy balance on the steam-side provides the total heat 


transfer in the boiling section, 


where h is the enthalpy of saturated vapor at the 


3 
pressure specified. The required average gas temperature 


entering the boiling section is 
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3. In the superheating section, the steam temperature is 
raised from the saturation temperature to the outlet 
temperature specified in the initial conditions. The heat 


transfer in the superheating sections is 


Osh 


where Reg is the enthalpy of steam corresponding to the 
steam outlet temperature Teg: The gas inlet temperature 
+ is known from the initial conditions. 


gl 
The "pinch point", for the purposes of this model, is 


defined as the difference between the gas temperature enter- 
ing the heating section and the water temperature leaving the 
heating section (saturation temperature). This temperature 
difference will normally correspond to the smallest tempera- 
ture difference between gas and fluid in the heat exchanger. 
The only other likely location for the smallest temperature 
difference between gas and fluid is at the gas inlet, or 
Superheater outlet. Figure 3 depicts the typical gas and 
fluid temperature distribution in the WHRU. 

The model must provide some internal control over the 
pinch point temperature difference as defined above. This 
is desirable because the gas-fluid temperature difference 
at the fluid saturation point is indicative of the amount of 
heat transfer area which will be necessary to heat the in- 


coming water to the saturation point. As the temperature 
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difference at the pinch point increases the heating portion 
of the heat exchanger becomes more efficient from the point 
of view of heat transfer area required. The minimum pinch 
point temperature difference is set at 25F. After the 
interim temperatures are established the pinch point temp- 
erature difference is checked. If the difference is less 
than 25 F, the gas temperature at the heating section inlet 


is reset to 


From the overall energy balance 


nit. — « 
pg g | gl ga) 


and the energy balance across the heating section 


Ch Oa. Gre. op = 
pg g | g3 g4) 


a revised gas outlet temperature may now be calculated from 


enicisees ol 
T 54 L =a 
where 
na: ee a 
Neg fl 


28 






— ween 
arn 


+ gntog doaitq eit bedelidasee eis geass 
I et eoretetiib eds QI .beasloetio ef 


Rita 





- ’ + = 
r a — 
as 
e 
~~ rr 7 ~~? - * | 
i the S TPrsas ile 
; om La a) ; 
3°" a ad = » 4) fi 
da + a » }* BD 


osz priczsed eit atozos sonsisd wens or 


a 





After this new gas inlet temperature is established, all 
initial calculations are performed again, to establish a 
new water/steam mass flow rate and to recheck the pinch 
point temperature difference. Once these calculations are 
complete, the geometric parameters for the heat exchanger 


must be specified. 


B. GEOMETRY 

A fin-tube with helically curved extended surfaces on 
circular tubes is assumed for this model. The fins are 
segmented. The tubes are configured in banks of one row 
each and the rows are staggered. The base tube surface is 
taken from Ref. 7. The description of the finned tubes 


is as follows. 


dq. = tube inside diameter = 1.86 in. 

dq, = tube outside diameter = 2.00 in. 

dq. = fin root diameter = 2.00 in. 

Ne = fins per inch = 5.94 

Q = fin height = 1.015. in. 

oes = wength On: Cuc from £in tip. = 0.82 in. 

d = fin outside diameter = 4.03 in. 

Saha, Fin thickness’ .=, 0.048. in. 

We Oy fin segment width = 0.17 in. 

Ng = number of segments in 360 degrees = 38 
The finned tube configuration is shown in figure 4. The 


tube length, L, and number of tubes per row N are 


CAL 
chosen by the designer according to the frontal area, A,, 
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Figure 4: Segmented Fin-Tube Configuration 


desired. The tube layout is shown below in figure 5. 

The center-to-center tube spacing in the transverse direc- 
tion is 4.50 inches. The equilateral, staggered tube 
arrangement used in this model leads to a spacing normal 
to the gas flow, Sat of 4.50 inches and a spacing parallel 
to the gas flow, 5. of 3.90 inches as shown in figure 5. 
The heat exchanger height is established by the number 

Of WHRU passes (rows) and the tube spacing parallel to the 


gas flow ut 


30 


a 


héviu Lidssse 2] 2Hpled vepatioze js0c 
o-losd& scut off Gan (ew0t) eseseg URW 





wi-ai5 oOssnemps2 <:4 saps? 


"WORs gi guoyvel sdos si? .beul 


ial 






Th23S538 Stl? selttao-o2-3esa080D en 


fstésoiivee ga? .eefsni O2.% 2f at 


cl lebom 2Lay at beer teemenaes 


% 
eafonl OF.6 Se , S .well wap eels os 









& 
+ 


n 
iso) Of.£ Se oe wool? e@p andZ oz 


' 7? i 
tt 


Figure 5: Model Tube Layout 


In order to establish the minimum gas flow cross- 
sectional area the total "blocked" frontal area, ADs of 


the heat exchanger must be calculated from 


Ay = Neji Se 


and the minimum gas flow area is 


The total inside area available for heat transfer per pass 
is 
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The gas side surface area available per pass for heat 
transfer is the sum of the fin surface area and the bare 
tube surface area per pass. The fin surface area per tube 


is calculated from 


Agin = Ne LIN. (2 he wot2 te. ae We te) 
i a2 2 
= x (dey ~ do) ] 
where dep = d = 2h: The bare tube area per tube is 
Abt = T dq LIF ai is dq t. Ny L 


Therefore, the total outside area per pass available for 


heat transfer is 


Asp > i/o Bein ee 


The cross-sectional area for fluid flow is calculated from 


2 
cL ee 


= pul 
nab a 4 
With the mass flow rates, terminal temperatures and heat 
exchanger geometry established, the remainder of the model 
may be solved for number of passes, actual interim tempera- 
tures, location of the fluid phase changes, and the gas 


Side pressure drop. 
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C. HEATING SECTION 

The gas-side Reynolds number for the heating section is 
calculated initially using the gas bulk temperature to find 
the gas properties. With this Reynolds number, a j-factor 
arc a tLriction factor, £, are obtained from polynomials fit 
mo the data for tube layout number 5 in Ref. 7 (fig. 5). 
The friction factor will be used later to calculate the 
gas side pressure drop in the heating section. The j-factor 
is related to the heat transfer coefficient se by the 


following relationship. 





z 2/3 
5 eae ee 
By introducing 
Ny 
S — 
€ Re E - 


N N 
= u p2/3 ¥ on : 
R Pp e y, ‘ 
eg, rg Rag Pig 
and 
= 7 Bo. Pee, 
u rg 
where 
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Therefore, a relationship may be written for the heat 
transfer coefficient as follows: 
K 
eg rg 
fe) 
The water-side heat transfer coefficient is calculated by 


using the Dittus-Boelter correlation 


where all properties are obtained at the bulk temperature 

of the water in the heating section. Using the tube wall 
resistance together with he and Aer the overall heat transfer 
coefficient for the heating section may be written in terms 


of the inside area as 





where 
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The fin efficiency, Nes is calculated from the expression 


es tanh ML 
Ne ML , 
where 
M = Yh P/KA 


and, if the fin is approximated by a set of rectangular 
d.--d 

strips extending from the tube wall, L = =. Now, 

the cross-sectional area, A, of a rectangular strip may 


be written as 


So, uSing known geometric parameters and the thermal conduc- 
tivity of the fin metal, the parameter ML may be expressed 
as 


ME = CC vh_ 
g 


where 
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Now, 


tcanh-C-sh-— Asn 
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2 vay op 


The number of passes required inthe heating section is 
calculated, using the effectiveness-NTU method, in the 
following way. 

1. An average pass effectiveness for the heating section is 
calculated from the following expression for cross-flow 
effectiveness with both fluids unmixed obtained from 


Ref. 8, 


ED = iL — exp [S42 -NCo) -1) 


where 


min pe £ 
Cc = ' 
@ e 
max pg g 
ol aH 
N = NTU = = ‘ 
min 
| at ieee 


Water and gas properties are taken at bulk temperatures. 
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2. An overall heating section effectiveness is calculated 


from 


eee Cys ee fC 
( p min max) n eer 
Lo =. 
foa 5 ‘ee Oe re ee 
( p min max) n ~ Quin 
1 oe € 
p max 


where n = number of passes. This expression was obtained 
from Ref. 9. The overall effectiveness formulation was 
actually derived with the condition that fluids are mixed 
between passes. It can be shown [Ref. 9], however, that 
the error is not large when the expression is used for the 
case where fluids are unmixed between passes. 

3. The expression for See is now used in an iterative 
fashion starting with n = number of passes = 1. Each time 
eas is solved, the gas temperature into the heating section, 


teenizs Eound from 
g3 


€ an oneg. L. =o fy 
S Oa Tman, ool max g4 
WM — 
g3 € eargeampemin © 
oa “min max 


2 a So ie Tha) 
a Cee (perss Py yy = 
oe Cain ae ai 


4. Once the inlet gas temperature is established for a 


particular number of passes, an energy balance may be 
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performed to solve for the total heat transfer which would 
take place in the heating section for that number of passes. 


( ) 


m- € 7 ane 
Mh g “pg g3 g4 
5. The enthalpy of the water at the outlet of the heating 


section for a particular iteration may be found from 


and this yields the temperature of the water at the outlet, 


r for the number of passes under consideration. 


EZ 
6. The outlet water temperature, Teor may be checked to 

see if it exceeds the saturation temperature. 

7. This set of calculations is performed until Teo exceeds 
the saturation temperature. This naturally means that 
boiling begins in the last pass of the current total number 
Se passes. Initiation of boiling will be treated in the 
boiling section of the model. Therefore, the heating section 
ends with the pass prior to that which initiates boiling. 

The results of the above calculations are: (a) The total 
number of passes contained in the heating section, (b) the 
new gas temperature at the inlet to the heating section TY3! 
and (c) the new water outlet temperature Te5° 
With these new heating section terminal temperatures an 


average outside wall temperature for the section may be 


calculated from 
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where ALj/ALS are total inside/outside areas for the section 
and Top/Tep are the new gas/fluid bulk temperatures for the 
section. The expression for To is derived from the follow- 


ing formulation for heat transfer in the heating section, 





Meiwep tS fh _ gb “wo 
at - R 
LRen ? 
which reduces to 
Ds 
Dae = ek = ob Eb? 
= th 
where 
i 
R = 
fe) Ne oe Alo 
and 
ib 
heat = , 
th ae: EL 
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This gas-side film temperature is now introduced at the 
beginning of the calculations for the heating section, to 
replace the gas bulk temperature. All calculations are 
performed again with Tf: 

The gas-side pressure drop in the heating section may 
now be calculated using the previously obtained friction 
Paector, £, from the correlation contained in Ref. 7. The 


gas-Side pressure drop is 


AP om cE Gnax Ne/e Hw) 0.14 
Q Uy 
where 

m 
G — 2bign 3 
max Be: 

min 
yee = number of tubes per row 
fe) = gas density 
We = gas viscosity at wall temperature 
Uy = gaS viscosity at bulk temperature. 


D. BOILING SECTION 

The boiling section is solved pass-by-pass since most 
of the section will involve two-phase flow on the cold side 
and the inside heat transfer coefficient will change with 
Gidlicy, x, and heat flux, gq". The correlation for the 


two-phase region inside heat transfer coefficient selected 
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for this model is that recommended by Tong [Ref. 10] for 
both nucleate boiling and forced convection. The equation 
for the two-phase flow heat transfer coefficient was 


given by Schrock and Grossman in Tong [Ref. 10] as 





h 
" n 
iy, fg 7G 
The constants are given by Wright [Ref. 10]: B= S70 win 


A = 3.5 ane age n = 0.66. The heat transfer coefficient 


assuming a total liquid flow is 


K d. G(1l - x) C u 
Ao = 0.023 a1 = Sa 
i Uo Q 


and the Martinelli parameter is 





The limits of the data for this correlation are 


quctatys, »O.05— 0.57 


heat. EDusx: go son0- = TAS TO” BTU/hr-£t*. 


In this model, heat flux in the boiling section will be at 
levels considerably below aie BTU/hr-£t- and the full 
range of quality must be considered. Therefore, in order 
to observe the performance of the Schrock and Grossman 


correlation at lower heat flux and quality above 0.57 a 


plot was made (fig. 6) of the ratio he pe/ Bp vs. quality 
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where hy is the heat transfer coefficient from the Dittus- 
Boelter correlation for water at saturated liquid condi- 
tions. The correlation of Dengler and Addoms as given in 
Tong [Ref. 10] for entirely forced convection vaporization 
tea ineLuded for comparison. Tong (Ref. 10] indicates that 
low heat flux and high water-vapor mixture velocity favor 
the forced convection mechanism. After a review of the 
information available in Tong, it is not entirely clear 
whether the forced convection mechanism or the nucleate 
boiling mechanism is dominant for the relatively low flow 
velocities and low heat flux of this model. In any case, 
the Schrock and Grossman correlation applied at low heat 
flux represents a type of compromise between the entirely 
forced convection assumptions of Dengler and Addoms and 
the "mixed" assumption of Schrock and Grossman. When the 
Schrock and Grossman correlation is applied under the 
Sonditions of very low heat flux and quality (x < 0.05), 
as in the first pass of the boiling section of this model, 
it can be seen that the ratio Ape he becomes less than 
one. This is both unrealistic and computationally diffi- 
cult to handle. Therefore, in the boiling section, the 
model assumes that Apt = hy for the first pass or for any 
pass where x < 0.05. This assumption will produce a more 
realistic but still conservative design. 

One other simplification should be explained prior to 
discussion of the actual model calculations for the boiling 


section. Since the boiling section calculations are performed 
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for one complete pass at a time, the average quality, a in 
a particular pass is used as the "local" quality to compute 
the inside heat transfer coefficient from the Schrock and 
Grossman correlation. 

The first pass in the boiling section, as indicated 
in the calculations for the heating section, will, in 
general, involve both heating water to the saturation tem- 
perature and the initiation of boiling. From the calcula- 
tions for the last pass of the heating section we have the 
following quantities for the first pass of the boiling 
section: water inlet temperature we, oy 
outlet temperature (pt ). The gas side heat transfer 


out 
coefficient, Seu is obtained, using the average gas tem- 


average gas 


perature in the boiling section, in the same way as for 
the heating section. 

A rough estimate for the gas temperature into the first 
pass must be obtained in order to begin calculations for 
eat pass. This is accomplished by calculating a rough 
overall heat transfer coefficient, dao! for the pass using 
the tube wall resistance, he? and neglecting inside resis- 
tance. A pass effectiveness, co is then calculated assuming 
that only boiling takes place in the pass. The expression 
for effectiveness for fluids unmixed between passes for 


C = 0 becomes 


€ = (1 - exp(-N) 
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where N = NTU. Using Ep calculated from this expression, 


the rough gas temperature for the first pass is 


Soe 1 


oi pee Tein ~ "Sout 
Jin e - 1 
p 


Since this pass has been simplified by allowing the inside 
heat transfer coefficient to be computed by the Dittus- 
Boelter correlation this initial guess for Tr Ls: not 
essential. It would, however, become ieeeeon if we were 
attempting to "split" the pass with two inside heat transfer 
coefficients, one for heating and one for boiling. Instead, 
the inside heat transfer coefficient, et is calculated for 
the entire pass uSing the Dittus-Boelter correlation with 


water conditions corresponding to the average water condi- 


Gieons. in the heating section of the pass, that lis, 


With this inside heat transfer coefficient, an expression 
for the overall heat transfer coefficient for the first 


pass may be written, 


Sethe iE 

Ow l x A; In(d, dq. Ai l 
baa oaniky N G A 1 
he w t/p op AG Ny 


The pass effectiveness, ene is now calculated in the same 


; ; it 
Manner aS previously described. With ED! the average gas 
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temperature in is calculated from 


I aL 4b 
ot Beemer ey mae Souk 
Jin I Cc - 


. Cc 
p “min max 


The pass heat transfer may be calculated from 


In order to determine the location of the interface between 
heating and boiling the total inside area devoted to heating 


in the first pass is calculated from 


nt 
al ,DUS, satf Fin 
mot ik ip 
h _ 
E 5 
out in 


Subsequent fluid passes in the boiling section are 
calculated until the quality at the outlet of a pass is 
found to be greater than 1.0. When this occurs, the boiling 
section is allowed to end with the last pass for which 


x < 1.0. A “mixed" pass where both boiling and superheating 


46 










"i caso Secs moza yoelsdtos, Jel suo sot soisses: 


a few fa 


=f =) ii a0 


most betsivules i yiilewp selsus 


- 
« 


aia 


i et o: a i 

: _ : zue 
‘sac 2587 ter: et Fo intjenol add eniersegeb oF 
tise o+ becoveb egné ehden? faves ofa patiiod bra p 4 


fi 
mez? Seseissies af eenq sett? | sci 


act ,eTue2c effio «ate «08 
is:dw oct aeog teak eda dear 
ae 

‘2 oh: pailso@ sad eaete 6s 


takes places will be calculated in the superheating section 
of the model. 

After the first pass, boiling section fluid passes are 
calculated in the following manner. 
1. The steam quality at the outlet of the previous pass 
becomes the inlet quality for the current pass, Xone The 
average gas temperature at the inlet to the previous pass 
becomes the outlet average gas temperature, T , Lor the 

Jout 

current pass. 
2. Using the last overall heat transfer coefficient calcu- 
lated for the previous pass, aa approximate values of 
NTU and ED are calculated. 


3. With this approximate ED the average gas temperature 


into the current pass is calculated from 


Go. ok = Gh 
T = PB Fat Jout 
Jin ase) =o 
4. Using this gas temperature in, T , the pass heat 


bial 
transfer, Qo is calculated. This allows the calculation 


Sm enthalpy out,. h and Gqualeey Oue, x as well as 


f i: Cue 


out 
the average quality in the current pass, a 
5. The inside heat transfer coefficient, Ape’ may now be 
calculated, 
QO /A. 
hog = Blea, + Be) oh hy 
P fg ei 
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where the Martinelli parameter ga is calculated using 
average quality, Xo! and all stkor terms are as previously 
described. 
6. A new overall heat transfer coefficient is now calcu- 
lated from which a new pass effectiveness can be calculated. 
7. This new pass effectiveness yields a new average gas 
temperature into the pass from which a new pass heat 
transfer may be calculated. 
8. The new pass heat transfer is compared with the pre- 
viously calculated pass heat transfer. If the two heat 
transfer calculations do not agree to within 5% the entire 
set of calculations is repeated. 
9. When the pass heat transfer calculation converges, the 
final Quality at the outlet of the pass. 1s calculated. If 
the quality is less than 1.0 another pass is calculated for 
the boiling section. 
ieee ce the quality is 1.0, the boiling section ends with 
the current pass. If the quality out exceeds 1.0 the boiling 
section ends with the previous pass. 

Calculations for the average gas film temperature and 
gas-side pressure drop are performed in the same manner as 


for the heating section. 


E. SUPERHEATING SECTION 
The first pass of the superheating section will, in 
general, involve both boiling and superheating. The inside 


heat transfer coefficients will be quite different for the 
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boiling and superheating portions for the pass, and, for 
this reason, the pass will be "split" to allow more accurate 
calculation of the heat transfer in these two regions. An 
outline of the calculation scheme is as follows: 

1. Calculate the gas-side heat transfer coefficient for 

the section in the same manner as for the previous two 
sections. 

2. Estimate the average gas temperature into the first 
pass. This is only a rough guess for Le , to begin the 
calculation. = 


3. Calculate the required heat transfer in the boiling 


portion 


where he is the enthalpy out of the last pass of the 
in 


boiling section and h is the enthalpy of saturated 


satv 
vapor. 

#. Using Ob and the average quality in the boiling portion 
of the first pass, calculate the inside heat transfer 


coefficient, h using the Schrock and Grossman correla- 


EDL 
tion, as in the boiling section of the model. This, together 
with the gas-side heat transfer coefficient and the wall 
resistance, will yield an overall heat transfer coefficient 
for the boiling section, Dee: 
5. The calculation of the area required for the boiling 


portion of the pass is performed in an iterative scheme as follows. 
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a. An initial estimate for the required boiling area, 
Al, is made. 
b. The mass flow rate of the gas over the boiling 


portion is then 


c. An effectiveness for the boiling portion may be 


calculated from 


= £ rb 
b b 
cata (8 4 = ) 
wen Sin sat 
where 
b aS 
ae = m pg” 


d. Using the previously calculated overall heat 
transfer coefficient, the estimated boiling area and the 


consequent gas-side heat capacity, NTU is calculated from 


e. A "test" effectiveness, 2 a can be found by applying 


this NTU in the relationship for effectiveness 


Ce es l = exp(NTU). 


' 3 : 1B 
Es The two effectiveness calculations, Eb and En! are 


compared. 
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CL) Sis ee is greater than Eye then the original 
estimate for the boiling area, Al, was too large for the 
heat transfer required and the specified gas temperature 
into the pass. So, the estimate for the boiling area is 
decreased and a new Al is established. 

C20 are ee is less than Er Al is increased. 


g- A new gas-side heat capacity is calculated from 


ch = = mc 
ip G Pg 


h. A new boiling portion effectiveness is calculated 


using the latest boiling area estimate 


in sat 

i. As before, an NTU is calculated with the new Al. 
This generates a new "test" effectiveness, ep which is 
compared with Ene EE ey and es are not equal, a new Al is 
calculated based on the comparison. 

j- These calculations are continued until oa = e,- 
The current Al is then set equal to Aen 
6. Once the area required for boiling, cee for a particular 
average gas inlet temperature has been established, the 


following first pass quantities can be calculated: 


aie cA. = area available for superheating 


sh ne b 
3 ip 


aL 
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b. mt = gas flow rate over the boiling portion 
(og in” = gas flow rate over the superheating portion 
d. 7p? = gas temperature out of the boiling portion 
Jout 
7? = 7 - weet 
Sout Jin me € 
Gg Pg 
e. sae = gas temperature out of the superheating 
Jout 
portion 
A; Ah c= AL 7? 
psh s P Sout Jout 
= 
Jout aot 
x 
, sh b 
(The weighted average of T and T must equal 
out out 
the T previously calculated.) 


out 
7. Since the heat transfer to the cold side must satisfy 


an energy balance involving the average temperatures of the 
gas, the expected heat transfer in the superheating portion 


may be calculated from the heat release on the gas side 


Q = mE. €. (Cr aa i 
p ep ie SI Sac 


and the required heat transfer in the superheating portion 


as 
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8. The enthalpy and temperature of the steam out of the 
first pass may be calculated from the assumed heat transfer 
in the superheating portion, Qon° 
9. The steam properties for the superheating portion may 
now be found. These properties, along with a and wall 
resistance, will yield an overall heat transfer coefficient, 
aaa for the superheating portion. The Dittus-Boelter 
correlation is applied to calculate the inside heat trans- 
fer coefficient. 

10. The heat capacities on the gas and fluid sides are 
calculated, using the current gas mass flow rate over the 


superheating portion, fie to calculate the gas heat capacity. 


The heat capacities along with wag and ach yield the effec- 
tiveness for the superheating portion of the pass, E Sh’ 

ll. A new average gas temperature into the pass may now 

be calculated from the superheating portion effectiveness 
and the superheating portion heat transfer required to 


satisfy the original heat balance from the expression 


Q 
a a sh 
: in - ie ‘ € on 
sh “min 
. * . 
This temperature, T , 1S a measure of the temperature 
in 


which would be required to produce the originally specified 
Ooh under the actual conditions of the heat transfer char- 


acteristics of the superheating portion. 


38 
















vie *o too mkede edt Se siptersgies bas yylad 
‘OF Smeas cen bemveds Sit mom’ petelesign ed 
-4g0 »02920¢ ese 


we » 


‘‘osiek 


on awLe2og teistass2aqua sat rot esittscory, 
: : 
leu Drea: ai ftiw gnols ,eeidasgott seefT «<dnvet od 
(Iso) tethers 3sed. Iiszeve as blety fitw ,onns: 
JisOd-2u37id. gf  -.aeigeeq. patseatissage edd ‘ 
Pets 32 silent s&¢ sadeleolep of bsilage 1) noksaia 
. 2notositeos 
4bie Biol? Bre ger ont 20 2e 
SI WOlLl #88m aap tn9tmeo Ss prem ber siog 
Sed esr eis &9atualeo os Pat ,AOLO 308 OA. eri 
Tt : te : 
etis eff Blaby “7Asboe pes iasy poole (BS 2538gEo 


7 .eesg ads ico noiseaq prigebiiague eae San 


ac eteo sdy otf expyaveumed. 6b Sestere wen ee 


1 2 


S2SOh°VLISCSESS BOI IOg DISS eiieegie Cit moss Beseiow 
a ‘igpes TSieGgets Jaoh NOLstog Peisemepegue 


J#2etqxs shits m2? sodsiad shee ideteion ails 


- & 
0 SInesem 5 SI~, T 


st De = "2 z 
sce yllsnspruc. sd) soybome es Bea ac biwow ae 


_ eT? Jeec sit to enpbsibnos wlevsoe ee Sh 


em = 





* 
we Lt the T is the same as the originally specified 
alli 
average gas temperature into the pass, T , then the heat 
Jin 
sh. 


under the current distribution of area between boiling and 


transfer in the superheating portion, Q is possible 
superheating and the heat balance is satisfied between steam 
and gas sides for the specified gas temperature in. [In 
this event, the calculations for the first pass of the 
superheating section are complete. 

oe ILE tT is less/greater than the current T _ the 

in in 

average gas temperature into the pass is increased/decreased 
and the entire set of calculations is performed again until 
* sé 
Baio in 

14. The second and subsequent superheating section passes 
are calculated in much the same way as in the heating section. 
That is, an overall pass effectiveness is calculated for a 
specified number of passes, using the multipass effectiveness 
relationship from Kays and London [Ref. 9]. 

15. The average steam-side heat transfer coefficient, Des 
is calculated using the Dittus-Boelter correlation with 
steam properties at steam bulk temperature for the remainder 
of the superheating section. The overall heat transfer 
coefficient and the average pass effectiveness, ete are then 
calculated. 


16. The formulation for overall effectiveness with the 


number of passes, n, variable, is 
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As in the heating section this expression is used in an 
iterative fashion, beginning with n = l. 
i7.. Each. time By is solved, the gas temperature into the 


superheating section is calculated from 


where Tr is the temperature of the steam into the second 
in 
pass of the superheater and T is the average gas tem- 
out 


perature out of the second pass of the superheater. This 
gas temperature in is compared with the originally specified 


gas temperature into the heat exchanger, Eas 
18. If the calculated gas temperature in, T ats is less 

than eral the number of passes is increased aces and the 
calculation 1s performed. .again« If T exceeds Teer 
superheater is allowed to stop with he eae pass. 


the 


Calculations for the average gas film temperature and 
the gas-side pressure drop are performed in the same manner 


as for the heating and boiling sections. 


F. MATCHING SUPERHEATER GAS INLET TEMPERATURE 
The model calculations thus far lead to a waste heat 


recovery unit design for which the gas temperature at the 
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inlet is, in general, lower than that specified in the initial 
conditions. This condition occurs because the superheater 
ends with the last complete pass calculated for which the 
entering gas temperature is less than or equal to the 
originally specified WHRU inlet temperature. This outcome 

is, of course, unsatisfactory since the inlet gas temperature 
is an independent quantity which is a function of the horse- 
power setting on the gas turbine. An iterative scheme is 
employed for matching the T calculated in the superheating 


in 


section with the T specified in the initial conditions. 


gl 

Since the model does not allow the calculation of frac- 
tional passes and since the scheme for the calculation of 
superheater passes will not allow inclusion of a pass which 
would require a gas inlet temperature higher than that 
specified in the initial conditions, the final gas inlet 
temperature calculated for the WHRU will generally be lower 
Bian that initially specified. In order to match this 
Ssuperheater gas inlet temperature with that specified 
initially, the heat exchanger gas outlet temperature is 
allowed to rise. Each time the gas outlet temperature is 
increased, the entire waste heat recovery unit calculation 
is performed again, and the calculated superheater gas 
inlet temperature is compared with the gas turbine exhaust 


temperature. When a match is achieved the design is 


fixed. 
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G. OFF-DESIGN CALCULATIONS 

Once a particular design has been selected the designer 
must test that design at several off-design points. A 
design will define the following WHRU characteristics: 

Dimensions (length, width, height) 
Operating pressure 

Superheater outlet temperature 
Fin-tube configuration. 

The off-design calculation requires that the selected 
WHRU design performance be predicted for a gas inlet tempera- 
ture and flow rate different from that which was used to 
produce the selected design. The design model can be used 
in an iterative fashion to make this calculation. All 
model calculations are made as prescribed in the previous 
sections, with the exception of the pinch point calculations 
which are not performed. 

The off-design procedure requires that the designer 
fix the WHRU frontal dimensions, operating pressure, fin- 
tube configuration and minimum gas outlet temperature. A 
design is then produced for a gas inlet temperature and 
flow rate at an off-design point. This design is checked 
to determine whether it conforms to the height dimension 
of the WHRU designed at the gas conditions of the off-design 
point. If the height dimension is not matched then the 
number of WHRU passes must be adjusted until a match with 


the design-point WHRU is achieved. The number of WHRU passes 
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may be altered by either adjusting the WHRU gas outlet 
temperature or the superheater outlet steam temperature. 

The effect of either of these adjustments is to change 

the steam flow rate and the heat transfer rate. It would 

be possible to adjust the steam flow rate and allow the 
Superheater steam outlet temperature to remain fixed. 
However, for the method of this model, the computation time 
required for this means of control would be large. The 
actual procedure used produces designs relatively quickly 
which closely approximate the results which would be 
obtained by a fine control on steam flow. Each time the 
Superheater steam outlet or gas outlet temperature is adjusted, 
a new design is produced. This procedure is continued 

until a design is produced for which the height (number of 
passes) equals that of the design point WHRU. When the WHRU 
height at the off-design gas conditions is equal to that 

of the design point WHRU, the gas inlet temperature is 
matched with the gas inlet temperature specified for the 


off-design point, using the procedure described in Section E. 


H. COGAS SYSTEM OUTPUT MODEL 

The COGAS system output model calculates performance 
parameters for the combined steam and gas turbine system for 
a Specified gas turbine input power. The WHRU model pro- 
vides a steam flow rate, pressure, and temperature for a 
specified set of gas turbine exhaust gas conditions which 


correspond to a particular gas turbine horsepower setting. 


58 


a » ad 
. 
a = rs Se 
=a) 
+ 


in 


* 


© 
— 


=j 


inv! ats 4 oisnatbn rand f= ve die 


: 
‘gs 4o£3 wo Feassadisg 


ha] 


a 







y bie, ‘ 


ws 


Ry 


a a fas oT ae 


= —-" 


sufbseeei2 to yerndis 


_ 


> teen od¢ Bos esses wold a 
I3 meerée sty Seve 
SoS Bisa 1 SsS.c0059 We 
Of efna 40 bansent 
s3me> te Beam Bele 
smh Senyhouc baa « 
Iie ia? sJemizotgg 
leadiee anit 


c 
j 
a 
i 


i 
ey 
in 


- 
} - ne se 
“53 ad et a Souk 
- . 
cm ~ (aay =P +ewe 
Lg - = =" see Oe 
e 
, > ‘ 
: ab | i 
. 


an. 
4 a) T 
ype le 
= oh 
_ - 
< =“ 3<¢ 
me 
er a. 
ww dad ast 


Ls | 
ry 


This WHRU output is applied in a simple Rankine cycle to 
calculate the steam-side power output. The gas-side 
pressure drop calculated in the WHRU model is used to 
arrive at a revised gas turbine horsepower. By combining 
these two power outputs, a COGAS system power can be calcu- 
lated. Other performance indicators such as specific fuel 
consumption, thermal efficiency and steam turbine share of 
the load are also calculated in the COGAS system output 


model. 


me) Rankine Cycle 


The following conditions are assumed for the Rankine 
cycle calculations: 
a. No line losses 
b. Steam turbine efficiency = meea= 0.85 
c. Condenser pressure = 2.0 psia (4 in Hg) 
d. Feedwater heater pressure = 15 psia 
e. Feedwater temperature = 200 F 
f. Pumping power required is neglected 
g. Fuel LHV = 18400 Btu/lbm 
The following input conditions are received from the 
WHRU model: steam turbine inlet pressure (Pi), steam turbine 


inlet temperature (T steam flow rate (me) gas turbine 


le 


exhaust pressure drop ease From Py and Ty the turbine 


inlet enthalpy (h,) and inlet entropy (sj) can be found. 


al 


The turbine exhaust steam quality, assuming isentropic 


expansion, x is calculated from 
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where 


S¢ entropy of saturated water 


S54 entropy increment for evaporation. 


The steam turbine exhaust enthalpy assuming isentropic 


expansion, Ao: May now be calculated from 


where 


> 
ll 


enthalpy of saturated water 


> 
i} 


fg enthalpy increment for evaporation. 


The isentropic turbine work is 
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In this model it is assumed that heating steam for the 
feed-water heater is provided, via a reducing station, from 
the main steam line. The fraction of the steam mass flow 


rate from the WHRU required for the feedwater heater is 


h -h 
a fw 2 
hy ~ Ay 


where 


a 
iH] 


enthalpy of the feedwater 


2) 
i} 


enthalpy of the condensate at feedwater 
heater pressure. 


Therefore, the steam turbine power output is 


2. Gas Turbine Performance Calculations 
The original gas turbine horsepower, exhaust gas 
flow rate, and exhaust gas temperature inputs to the WHRU 
model were hand-calculated using figure 4 of Ref. 11 and 
figures 1, 6, and 8 of Ref. 10. The assumptions used to 


calculate the WHRU model inputs were: 
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Ambient temperature = 100 F 
Ambient pressure = 14.696 psia 
Humadtey, =  0..0 


Fuel LHV 


18400 Btu/lbm 


inter toss = 4.0 in. HO 


Exhaust loss = 6.0 in. H,0 
In the COGAS system output model, the gas turbine horse- 
power correction factor for the gas-side pressure drop in 
the WHRU (BHP for 6" loss/BHP for WHRU loss) is found by 
solving 


Chhp = 21-0125 + 0.002125 APL. , 


where Chop = BHP for 6" loss/BHP for WHRU loss. This rela- 


tionship was found by extending the duct loss relationship 
Mmuperg. 7, Retr. Li linearly from the point of 6 in. HO 
exhaust duct loss already assumed. The final gas turbine 


horsepower is 





* 
where ae is the input gas turbine horsepower. The total 


COGAS system horsepower, P is found by adding the steam 


Tom 
turbine and final gas turbine horsepowers. 
In order to calculate the gas turbine specific fuel 


consumption fig. 4 of Ref. 11 is entered with the final 
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ca 


gas turbine horsepower, and the s.f.c. is read from the 
propeller law curve superimposed on that figure. The s.f.c. 
is 


correction factor for exhaust loss in the WHRU, Core! 


found in the same manner as C Loom tig. 6 of Ret. 11 


bhp 
where 


Occ = (s.f.c. for total exhaust AP)/(s.f.c. for 6" exhaust AP) 


and the gas turbine specific fuel consumption is 


The specific fuel consumption for the combined system is 
calculated by first finding the fuel consumption rate of 


the gas turbine 


oe ae 


Meyel 


mme COGAS system s.f.c. is then calculated from 


a i Meuel 
COGAS Bae 





To allow the direct comparison of the COGAS system with the 


all gas turbine system the specific fuel consumption of the 
ry 


gas turbine at the new, higher COGAS system power, Si Cwnet 


is calculated in the same manner as Sebi e: 
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The thermal efficiencies for the gas turbine at 
input power, the COGAS system, and the gas turbine at COGAS 


system power are calculated from the relationship 


es 2545 
Veh Saboc. LAV 


Additionally, the steam turbine share of the load 


is calculated from re are aa 
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III. RESULTS AND CONCLUSIONS 


A. BACKGROUND 

The foregoing model was applied in a computer simulation 
program written in FORTRAN IV for the IBM 360-67 computer. 

A listing of this program appears in Appendix A. Addition- 
ally a set of supporting programs was written for water, 
steam, and air properties. A listing of this set of pro- 
grams appears in Appendix B. 

In order to test the model and develop an understanding 
of the behavior of the WHRU, a set of designs were produced 
for a variety of initial conditions. These designs were 
then reviewed, and an analysis of the model behavior was 
performed. After completion of the model analysis, several 
of the designs were selected for observation at off-design 
conditions. The off-design simulations were produced for 
gas turbine input powers corresponding to ship speeds of 
9, 16, and 20 knots, which represents the entire range of 
gas turbine input powers considered in this thesis. Finally, 
two designs were selected for detailed testing at off-design 
conditions. These latter two designs were tested at gas 
turbine input horsepowers corresponding to speeds of 9 to 


20 knots, at one knot increments. 


B. DESIGN VARIABLES 
The following set of design variables was selected to 


produce the basic set of designs. 
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1. WHRU operating pressure: Three WHRU operating pressures 
of 400, 600, and 800 psia were selected for design production. 
2. WHRU geometric scale: In order to test the model 

response to a range of fin-tube sizes, designs were produced 
for the basic fin-tube configuration described in the geometry 
subsection at 1.0, 0.75, and 0.50 scale. 

3. Frontal dimensions: The approximate WHRU space constraints 
for a DD-963-type engineroom are shown in figures 7, 8 and 

9. Consistent with these constraints, frontal dimensions 
meets x 12" and 12'" x -r5S"-were selected for the initial 
design set. 

4. Gas turbine input power: The gas turbine input horse- 
power determines the input gas conditions, temperature and 
flow rate, for the WHRU. Since the most advantageous use 

of the COGAS system was assumed to be for "cruise" speed 
conditions, the range of gas turbine input horsepower con- 
Sidered corresponds to the gas turbine power required for 

a ship speed range of 9 to 20 knots, with the gas turbine 
Operating alone. The three gas turbine input powers selected 


for the production of WHRU designs are shown in Table I. 


Gas Turbine Input Powers 


GTHP Approx. Speed (kts) Temperature (F) Flow Rate (lbm/hr) 
16421 20 849 407589 
8526 16 742 328641 
1684 9 689 £597 3d. 
Table I 
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5. WHRU gas outlet temperature: The WHRU gas outlet 
temperature is fixed at 400 F. Reference 3 states that 

the principal cause of "hot" corrosion in the WHRU would 

be due to the formation of H5SO, from an expected concentra- 
tion of S03 of up to 50 parts per million. Reference 13 
states that dew point condensation of vapor in the exhaust 
gas becomes a problem for corrosion if the stack temperature 
is depressed below 250 F. Therefore, a minimum average gas 
outlet temperature of 400 F is considered to be reasonably 
safe from the aspect of gas-side corrosion prevention. It 
should be further noted that the outlet gas temperature 
calculated by the program is an average temperature, based 
on the heat release from the gas. Therefore, the actual 
temperature of the gas in the stack will be cooler than 

400 F in some locations due to the non-uniform distribution 
of the gas temperature as it leaves the WHRU. It should 
further be noted that, although this temperature is a fixed 
design variable for design production, it will vary somewhat 
from 400 F in response to the program calculations for 
setting the minimum pinch point and for matching the gas 
inlet temperature with that specified in the initial 
Gonditions. 

6. Superheater outlet steam temperature: The target WHRU 
superheater outlet steam temperature is 650 F. This tempera- 
ture will also vary somewhat in response to steam flow rate 


changes brought about by the program calculations to match 
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the gas inlet temperature with that specified in the 
initial conditions. 
7. WHRU water inlet temperature: The WHRU water inlet 
temperature is fixed at 200 F. 
8. The other major COGAS system parameters are fixed as 
follows: 

Condenser pressure: 4 in. Hg 

Steam turbine efficiency: 0.85 

Feedwater heater pressure: 15.0 psia 


Fuel LHV: 18400 Btu/lbm. 


Ce THE DESIGN SET 

A set of 54 designs was produced for the combination of 
design variables just discussed. The array of design 
variable combinations is described in Table II. The summary 
computer output pages for the results of each design run 
are presented in Appendix C. 

As discussed in Section A, the purposes for producing 
this design set were to test the model and to promote an 
understanding of the likely behavior of a WHRU for a range 
Bae conditions. To this end, the following set of output 
variables was considered. 

1. WHRU output 

a. height (number of passes) 

b. gas inlet temperature 

€. gas outlet temperature 


,d. steam flow rate 
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e. gas side pressure drop 
f. pinch point temperature difference 
g- heat transfer rate 
2. COGAS system output 
a. gas turbine horsepower (including loss from the 
WHRU gas-side pressure drop) 
b. steam turbine horsepower 
c. COGAS system horsepower 
d. steam turbine share of the load 
e. specific fuel consumption 
GEL ivs shac sot eCOGAS system 
(2) s.f.c. of gas turbine at COGAS system horsepower 
f. thermal efficiency 
Chey om 


C2)rean 


th of COGAS system 


e, of gas turbine at COGAS system horsepower 
Clearly, several of these variables are related, such as 

aie. Gio and Ten? and no additional information is provided by 
studying all variables in a closely related group. Therefore, 
the following subset of output variables was selected for 
analysis. 

1. WHRU height 

2. Steam flow rate 

3. Gas-side pressure drop 

4. Pinch point temperature difference 


5. WHRU heat transfer rate 


6. WHRU gas outlet temperature 
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7. Steam turbine horsepower 
8. COGAS system horsepower 
9. COGAS specific fuel consumption 

The above output variables are presented in tabular form 
using the format of Table II. Table II also serves as an 
index for the set of basic designs. Additionally, the 
significant trends in the variables are presented graph- 
ically. Where the trend of the results was ascertained to 
be the same for the two frontal dimension sets, only the 
values for the 12' x 12' front were presented graphically. 
Where clarity was better served, only one scale was plotted 
as representative of the trend described. 
Berect Of Scaling. 

The only significant trends noted with respect to scaling 
were a decrease in WHRU height with a decrease in scale 
and the same trend for gas side pressure drop (see Tables 
Pheer and iV). The trends. are shown in figures 10 and ll. 

The basic fin-tube dimensions of the model were scaled 
by Eactors of 0.75 and 0.50. For any scale with constant 
frontal dimensions, the total inside and outside areas 
remain the same. That is, for a constant length tube, the 
scaled inside/outside area (AZ /A2) is proportional to the 
full scale areas (A, /A,) with the scale factor, s, as the 


constant of proportionality. Hence 
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Since constant frontal dimensions are to be maintained, 
the number of tubes per row is scaled up by the scale 


Bactor, so 


N 
s = t/p s £7 
; = A. = w/P 
Ain : SA, and op a A, , 
where Ni/p is the number of tubes per pass. Therefore, the 


scaled outside and inside areas for a pass for constant 
frontal dimensions are the same as the full scale areas. 
Now, the scaled cross-sectional area for fluid flow on the 


inside will be, 








which is equal to the Reynolds number for the unscaled 
geometry. It can also be shown that the outside Reynolds 
number remains constant with scaling. Recalling that the 
Dittus-Boelter correlation was used to calculate the inside 
heat transfer coefficient for the single phase regions and 
used with a multiplier for the two phase heat transfer 
coefficient and assuming constant fluid properties, the 


Nusselt number remains constant 
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So, the scaled inside heat transfer coefficient 





increases with a decreasing scale factor. Similarly, for 


the scaled gas-side heat transfer coefficient, 





an increase is seen with a decreasing scale factor. 
Therefore, in addition to a decrease in height attributable 
to a scaling down of the fin-tube dimensions, there is also 
an improvement in heat transfer which results in a further 
reduction in height by elimination of passes. This trend 
is supported by the results. 

Recalling the relationship used to calculate the gas- 


Side pressure drop, 


it can be shown that the scaled gas-side pressure drop is 


inversely proportional to the scale factor, 
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where C is a constant based on the assumption of constant 
gas properties between scales. For the same number of 
passes, then, an increase in gas-side pressure drop could 
be expected. However, since the number of passes is decreased 
with decreasing scale, a net decrease in gas-side pressure 
drop was experienced (Fig. 11). 

Some differences in pinch point with scale were observed 
in Table V, but no consistent trend could be established. 
It was expected that pinch point AT would remain fairly 
constant with scale for fixed pressure and input gas condi- 
tions. That result was apparent to some extent, particularly 
in the range of high gas temperature and low steam pressure. 
This finding is consistent with the manner in which the pinch 
point AT is established. That is, the program fixes an 
initial pinch point AT only when the rough pinch point AT 
feeless than 25 F in the initial calculations. At the lower 
pressures and higher gas temperatures this mechanism would 
not be operative. Also, the final and initial pinch point 
AT calculations are based on the average gas temperatures 
at the inlet to the pass where the AT is calculated. There- 
fore, the actual pinch point may be either lower or higher 
than the one calculated, depending on the actual gas temp- 
erature distribution at the inlet to the pass. The pinch 
point AT will be discussed at greater length in the analy- 


Sis of the effects of pressure. 
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Effect of Pressure. 

The most important effect of pressure on the model 
performance arises from the calculations which establish 
the minimum pinch point temperature difference. As related 
in the model description, the minimum pinch point AT is set 
at 25 F. The WHRU outlet gas temperature is adjusted in 
the initial model calculations in order to ensure this 
minimum AT. The initial model calculations establish the 
steam flow rate from a simple energy balance involving the 


gas inlet and outlet temperatures (T ), the gas flow 


gl’ "g2 


rate, the water inlet temperature ( ), the required steam 


ail 


outlet temperature (T and the steam pressure. Tentative 


e4) ’ 


interim gas temperatures for the heater inlet, T and 


gs” 


boiling section inlet, T are also established in this 


g2’ 


calculation set. In order to establish the minimum pinch 
point AT, the model tests the interim heater gas inlet 


temperature against the criterion 


where T is the temperature of saturated water at the 


£2 


pressure specified. If 753 is less than Te5 + 25 “then t33 


is set equal to T + 25 and a new gas outlet temperature 


te 
is calculated as follows, 
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As an aid to understanding the behavior of the model, it 

is useful to predict the WHRU gas inlet temperature at 

which this pinch point calculation becomes operative for 

a particular pressure. This is accomplished by establishing 
a "critical" heater gas inlet temperature, Ty for each 


pressure considered, 


where T. is the temperature of saturated water at the 


sat 
pressure specified. The ratio a can also be solved for 
each pressure considered, based on a constant water inlet 


temperature of 200 F and a constant superheater outlet 


temperature of 650 F, 


cei sient 
inet pnp, 


where Ae = enthalpy of saturated water at the pressure under 
consideration. Finally, assuming a minimum WHRU gas outlet 
temperature of 4-0 F the minimum WHRU gas inlet temperature 
at which the pinch point calculation comes into effect may 


be calculated as follows, 
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Table VI summarizes these calculations for each of the pres- 


sures considered in the design set. 


Minimum Gas Inlet Temperature for Pinch Point Calculation 


Pressure T. (F) O, T. (F) 
mipsialc J: eae rs A, cE 
400 470 s2L9 720 
600 497 ~262 770 
800 543 7299 878 

Table VI 


Now the pressure at which the pinch point calculation becomes 
Operative for each gas turbine input power considered is 
predicted. The pressure for which the pinch point calcula- 
tion begins to adjust gas outlet temperature is listed below 


for each gas turbine input power considered, in Table VII. 


Pressure At Which Pinch Point Calculation Takes Effect 


atte cower candi cc eee 
16421 849 800 878 
8526 742 600 770 
1684 689 400 720 
Table VII 
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It should be recalled that the outlet gas temperature is 

also adjusted by the iterative technique employed to match 
the model-calculated gas inlet temperature with the gas 

inlet temperature of the initial conditions, once the 

design is established. This is a relatively minor adjust- 
ment, however, and is never greater than 10 F. These adjust- 
ments to WHRU gas outlet temperature are apparent in the 
values of Table VIII and are displayed graphically for the 
@eeex 12" frontal dimentions and 0.75 scale in Fig. 12. 

The final calculated pinch point for each of the design 
sets reflected the influence of the initial pinch point 
calculation. These results were shown in Table V and Figure 
13. It is re-emphasized here that the calculated pinch 
point AT is useful only as an indicator of what the actual 
Pinch point might be for a particular heat exchanger design. 
This is true because the model neglects the possibility of 
a gas temperature distribution across the inlet of a WHRU 
pass and, instead, uses the average gas temperature for 
calculations. 

There are two other trends in the design results which 
emanated directly from the control exercised by the model 
over the pinch point AT and the resulting gas-fluid tempera- 
ture distributions in the designs produced. These trends 
were apparent when the effects of pressure on WHRU height 


and WHRU heat transfer/steam flow rate were examined. 
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The effects of pressure on WHRU heat transfer and steam 
flow rate are shown in Figures 14 and 15 respectively. 
Tables IX and X list the actual values obtained for WHRU 
heat transfer and steam flow rate for the designs produced. 
WHRU heat transfer and steam flow rate directly reflected 
the adjustment of the WHRU gas outlet temperature to main- 
tain the minimum pinch point AT. That is, as the WHRU 
outlet gas temperature was raised, the heat transfer rate 
and steam flow decreased. The small increase shown in 
Figures 14 and 15 in heat transfer and steam flow rate 
for 16421 gas turbine input power at the 0.50 scale reflects 
only a minor adjustment difference in T for matching 
Bc between 400 psia and 600 psia. cis 

The relationship of WHRU height to WHRU pressure is 
slightly more complex. Figures 16, 17, and 18 demonstrate 
the approximate gas-fluid temperature distribution along 
the total length of the WHRU for the 0.75 scale at the 
three pressures and input horsepowers considered. The 
general trend observed in these diagrams of decreasing gas- 
fluid temperature difference with increasing pressure is 
representative of all the designs produced. This decrease 
in gas-fluid temperature difference, taken by itself, would 
produce larger heat exchangers with increasing pressure. 
This trend may be observed for the 12' x 12' front at the 


/high gas turbine input horsepower in Figure 19. At the 


medium and low gas turbine horsepowers, however, the WHRU 
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height remained fairly constant or actually decreased with 
increasing pressure. This trend is explained by the 
decreasing heat transfer rate with increasing pressure 
attributable to the pinch point calculation adjustment to 
WHRU gas outlet temperature shown in Figure 12. 

The design set results for steam turbine horsepower 
output are given in Table XI. These results are presented 
Seaphically, for the 12" x 12° front in Figure 20. Two 
separate trends were identified. The steam turbine horse- 
power demonstrated an overall increase with increasing 
pressure for the high gas turbine input horsepower. Although 
there is a small net decrease in WHRU heat transfer rate 
with increasing pressure at the high gas turbine input 
power, the steam is rejected to the condenser at a lower 
quality for the higher pressures. Therefore, there is less 
heat rejection in the condenser for higher pressures. This 
better performance of the high pressure systems in the 
Rankine cycle more than compensates for the small decrease 
in WHRU output at these pressures. .For the medium and low 
gas turbine input powers the trend of decreasing WHRU output 
with increasing pressure (Fig. 20) was dominant in the 
Rankine cycle also, producing systems of decreasing steam 
turbine horsepower with increasing pressure. 

The design set results for the COGAS system combined 
horsepower are given in Table XII. The results are displayed 


Peepnicaiiy wa Figure: 21 tor the 12" x 12" front. It is 
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apparent that the system horsepower output trend followed 
that for steam turbine horsepower output. This was the 
expected result. This trend is modified only slightly by 
additional gas turbine horsepower losses due to a changing 
gas-side pressure drop across the WHRU. 

Finally, the system specific fuel Sonisameeien followed 
closely the results for steam turbine horsepower output. 
That is, slightly better systems were observed at lower 
pressures for the medium and low gas turbine input powers 
and a slightly better system was produced at higher pressures 
for the high gas turbine input horsepower. These trends 
are shown in Figure 22 and Table XIII. 

Effect of Gas Turbine Input Power. 

An increase in WHRU height with increasing gas turbine 
input power was noted in the results for the design set 
produced (see Fig. 23 and Table III). This increase in 
WHRU height generally followed the increase in WHRU heat 
transfer rate with increasing gas turbine input power, Fig. 
24. This trend is modified only by a slightly increasing 
outside heat transfer coefficient with increasing gas turbine 
input power (see calculations for gas-side heat transfer 
coefficient in the model description). 

As expected, the COGAS system specific fuel consumption 
improved (Fig. 25) with increasing gas turbine input power. 
This improvement represents the trend of improvement in 


the gas turbine itself with increasing power, since the 
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contribution of the steam system decreases with increasing 
power. 
Effect of Frontal Dimensions. 

The only significant trends for changing frontal dimen- 
sions were in WHRU height and volume (see Tables III and 
XIV). The expected result of decreased height with increased 
frontal area was observed. The gas-side pressure drop 
follows this decrease in height. The WHRU volume showed a 
Significant increase for increased frontal area for all 
designs produced. This result is consistent with the model 
formulation. That is, as the frontal area is increased 
and the tube length is held constant the area for fluid flow 
is increased with a consequent decrease in Reynolds number 
and inside heat transfer coefficient. The same result is 
seen in the outside heat transfer coefficient. As the 
frontal area is increased, the minimum flow area for the 
gas is increased, cauSing a drop in the gas side Reynolds 


number and heat transfer coefficient. 


Bs OFF-DESIGN RESULTS 

The off-design results were produced using the same 
computer simulation program used to produce the design set. 
The method used for off-design runs was described in Section 
F of the model description. As the computer program is 
presently designed, the designer must manually adjust the 
inputs to produce an off-design run. The procedure is 


reviewed briefly below. 
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First, the designer selects a WHRU design to be inves- 
tigated at off-design conditions. Once this design has 
been selected, the pressure, frontal dimensions, height and 
scale of the WHRU are fixed. The designer then enters the 
gas flow rate and gas inlet temperature corresponding to 
the off-design point to be investigated. A design is pro- 
duced with these gas conditions and the physical charac- 
teristics of the WHRU design under consideration. The 
designer then checks the number of passes (height) for the 
resulting design. If the number of passes does not match 
that of the selected design, the WHRU gas outlet temperature 
and/or superheater steam outlet temperature are adjusted 
to increase or decrease the number of passes until a match 
is achieved. Once the physical characteristics of the off- 
design point WHRU match those of the design-point WHRU, the 
gas inlet temperature is matched with the initial conditions, 
and the performance of the selected WHRU design is established 
for the off-design point under consideration. 

In order to test the feasibility of the procedure described 
above and to demonstrate the off-design performance of a 
wide variety of WHRU designs, ten WHRU designs were selected 
for off-design runs. These designs are listed in Table XVI. 
None of the full-scale designs was considered because the 
heights of the designs produced at that scale were, for the 
most part, either too large or borderline according to the 


Space limitation described in Figure 9. 
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The results of these off-design runs are given in 
Table XVII. The summary output pages for each of these 
off-design runs along with the design-point runs are pro- 
vided in Appendix D. In terms of specific fuel consumption, 
the performance difference between the design sets con- 
Sidered was not large. The performance range for fuel 
consumption is even less significant. Table XV demonstrates 


the performance range for the designs considered. 


WHRU Off-Design Performance Ranges 


fuel consump- approx. fuel consump- 
Power Ste ec’. tion (lbm/hr) tion difference (gal/hr) 
high BE 7AS od ei 0) 699250 — 7019.3 tort 
Medium .391 = .403 4429.2 = 4439.7 das 
low soo, — soc T764.6 = TIG7.5 0.4 

Table XV 


Since the fuel consumption comparison showed very little 
difference between the ten designs considered and since all 
designs considered satisfied the dimensional constraints 
of figures 7, 8, and 9, only two designs were selected for 
additional off-design runs. 

For these additional off-design runs, a pressure of 
600 psia was selected. This pressure did not produce the 
smallest designs but it is a steam pressure with which the 
Navy has some experience in terms of materials requirements, 


Maintenance, and operation. A scale of 0.75 was selected 
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which yields a tube outside diameter of 1.5 inches and a 

fin spacing of about 8 fins/inch. The rationale for this 
selection was that a scale of 0.50 with 1.0 inch OD tubes 
end a fin spacing of 12 fins/inch would probably be too 
susceptible to both inside and outside fouling. [In order 

to achieve the closest possible comparison, the design-point 
for both designs selected was for the gas conditions corres- 
ponding to the medium (8526 BHP) gas turbine input horse- 
power. Both frontal areas were used for these detailed 
off-design runs. Table XVIII summarizes the characteristics 


of the two designs selected. 


Final Design Characteristics 


design GT input 
run # design point scale front pressure height 
14 8526 Ge S P2312 600 3a2 
41 8526 O.75 ists 600 PaeM 


Table XVIII 


Table XIX summarizes the performance results for the 
COGAS system with the 12' x 12° front WHRU and Table XxX 
gives the results for the system with the 12' x 15' front 
WHRU. In both cases, the performance characteristics of 
the gas turbine alone, at the COGAS system horsepower, are 
provided for comparison. The summary output pages for each 


of these runs are provided in Appendix E. 
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The specific fuel consumption curves vs. brake horse- 
Meumes cor the 12° x 12" front and the 12* x 15' front are 
plotted in figures 26 and 27 respectively. The difference 
in s.f.c. performance between the two WHRU designs con- 
sidered was not significant. Geometric considerations were 
therefore used in selecting the better of these two designs. 
A geometric comparison of the two WHRU designs is given in 


Table XXI. 


Geometric Comparison of the Final Two Designs 


Frontal Total Out- 
Front Area Height Volume Side Area 
me x¢15' 181.9 Ze9 528.96 49586.4 
me x12" 144.8 Sre.2 464.64 2776.5 
Table XXI 


It was concluded that, since both designs met the height 
menstraint of figure 9, the 12" x 12" front was the 

better choice because of the significantly lower volume and 
total outside area requirements. 

After the selection of the 12' x 12" front WHRU as the 
final design, estimates were made of the possible COGAS fuel 
Savings for a DD-963-type destroyer over 1000 hours of 
Operation. The NAVSEC Standard Destroyer Profile (Fig. 28) 
was used to determine the number of hours a destroyer-type 


Ship spends at each speed during 1000 hours operation. 


1519 


at 
, 


-sazod sAs3d -.20 aveu2 nostmence Om > 
eta Saoxe "Or "SS add inp goss “Ste MS 


sonsxgitib sft .ylevizssceax To Boe OS 
“102 Sie. zeb UBHW ows ant neawsed 98 
iow, snalgaxebiéacm Siszenaee -tnsvebbepie 
-enplesb owt epeds Fé imsiee at eye ae 
ci gevio 2i enpiesh. UFEW ows end Jo noeigegueo 










< 
Ld 
= 
7 


eapieest ow: Lenid en Yo meaty sates ors: 


- 31/0 aye 
78ce2r oe. 82 e.$ 
saT7te boa igh mae 


smusoy tewol yisnsot?tapte eds Io sapeeed eolors 

étnetetiupes s6g~9 abt 00 
sit 2s OFHW snOa2 "SE ~ “Cl oho 2p webvonien ae 
2AD02 eldiseng sda fo shan etew eavamiten ,ft 


ae 
} x. 
~~ = 
to exyod 000! sevo .~zsyotteam SgeseEse-Oe 2 203 a 
esiizoz9 se fotreso p= slbauss2 DSB VAU edt ahah 
yy y7-Teyorteeb, B exved G6 uedaun ae scimzezeb ¢ bee 


von toh palaw& sine, FT) ne \ 


ee 


9c AYNODIA 


Galt x ein SHO 


ia] — 
Q @ 0} A Nw 9g Le 8] oO A tw Qg 


CLG ils ei oD cE SN sd a GHOO 


stp 


3 


CHH-dHe Wa 1) 


120 





4 





—~ o ee ee ———— oe 


of v “4 2 


~ — — —e 


38! KX GHH 2AA02 


ay, stp iu 





ae 


a 


Le duno 


£40 DM dng? Se00) 


N N hw w -— _ Lom m - 
Lap) A i) aS) ies) oO & nw g c6 6) o bh tw 9g 
eee eee — eee eee oe oe oT E 


CNOA) do iG hee? bl ee Os ai dei OU IO) 


ved 


<2. 


CAH-dHa/Wa 1) 


nieeBk 


S 


write 





WWWN Ww 
vue CO OD 
O° OOT GT 
G36 v*O 
T° 86 S‘*O 
9°46 4°0 
6 96 ae 8 
9 S6 G°T 
le 76 mn 
8°68 SZ, 
ARS 6°8 
vel AL 
i, SQ ood 
G* eg 9°9 
@ tS BQ 
a Sv a 9 
Oo Ge CL 
G TE O°ZL 
6° ve L 9 
ie aE 69 
& TT & TE 
LNAOWAd AW Id 


SATLVINWNASD LNAOWAd 


8c HUNOIA 


CSLONA ds, OAaeS 


OT-t 
SLONM 6 


ot eel el A Oe IS AT LAS ABN e 


AW NOTSSIW 30)JNa30dad 


122 


TUL 


yr 


+) 


bEMCEML Ok 








ca 


Figure 29 was used to obtain the shaft horsepower require- 
ments for each speed. Finally, figure 30 was used to obtain 
the s.f.c. of the gas turbine operating alone. 

Two estimates of fuel savings were made. The first 
estimate assumed that the COGAS system would be operated for 
speeds of 5 to 20 knots and would be secured outside that 
range. The second estimate employed the COGAS system for 
speeds of 5 to 23 knots. For the latter operating range, 
the criterion for the high COGAS speed was an imposed maxi- 
mum gas turbine horsepower input of 20000 BHP. Both fuel 
Savings estimates are based on the following assumptions 
and simplifications: 

1. The COGAS mode of operation implies the use of one 
engine (COGAS mode) on one shaft with the other shaft 
dragging. 

2. The gas turbine is assumed to operate at idle (1000 BHP) 
at speeds below 8 knots. 

3. Maneuvering combinations, when two main engines would 
normally be on the line, are not considered. 

&. In the speed range of 21 to 27 knots, in the pure gas 
turbine mode, one engine per shaft is used, with two shafts 
on the line. 

SB. From 28 to 32 knots four gas turbines are on the line. 

Table XXII provides the estimates of s.f.c. and fuel 
consumption for COGAS between 5 and 20 knots and for pure 


gas turbine for all speeds. The total estimated fuel savings 
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COGAS PERFORMANCE: 
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IDLE 
at 
10060 
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2205 
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4316 
5579 
1368 
8947 
10421 
12316 
14210 


16526 


3 
5 to 20 KNOT RANGE 


COGAS 
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TABLE XXII 


Page l 
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CONS 
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CONS 
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COGAS GT COGAS GT 


Speed Hours BHP SECS Set... CONS _ CONS 
a ihe. 17368 -519 676049 
22 43 20000 495 425700 
23 LS 22105 -482 159819 
24 r3 25263 - 463 152058 
25 7 28684 -448 ; $9953 
26 5 SZ6sl -435 T0972 
Zo 4 37895 -418 63360 
28 3 44210 - 482 62928 
7) 2 51580 - 460 71180 
30 3 58947 444 73517 
a1 3 69474 ~427 88996 
32 3 78947 ~415 98289 


4774897 5429976 
12.1% savings 
over 1000 hrs operation 
655088 lbm savings 
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over 1000 hours of operation was 12% which is roughly 
equivalent to 96000 gallons of distillate fuel. 

The estimates of s.f.c. and fuel consumption for pure 
gas turbine and for COGAS for the 5 to 23 knot speed range 
are given in Table XXIII. The total estimated 1000-hour 
fuel savings for this COGAS operating range was 19% which 
translates to approximately 152000 gallons of distillate 


fuel. 


E. CONCLUSIONS 

The waste heat recovery unit model and the accompanying 
Simple COGAS system output model provide the basic frame- 
work for additional studies of the application of the COGAS 
system to U. S. Navy gas turbine powered ships. The model, 
as presently formulated, provides a reasonable estimate of 
WHRU size required for the fin-tube configuration considered. 

Estimates of the COGAS system performance in the 5 to 
20 knot range and the 5 to 23 knot range are encouraging. 
Even considering the assumptions involved in making the 
fuel savings estimates, it is likely that a COGAS version 
of a DD-963-type ship would consume on the order of 300000 
to 450000 gallons less fuel in one year of operation than 
the current DD-963 class ships. 

The most probable COGAS system would be one which would 
be designed to operate in non-maneuvering situations for 
either speed range considered above. That is, a system which 


would be employed for maximum fuel savings in the cruise mode. 
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COGAS PERFORMANCE: 
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TABLE XXIII 


Page 
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COGAS GT COGAS GT 


Speed Hours BHP S.£.¢. s.f.c. CONS _ CONS 
21 “us 18421 - 340 — ye BA, 469735 676049 
22 43 20842 £335 495 300229 425700 
23 tS 23368 -324 - 482 113566 159819 
24 13 252263 - 463 £52058 
25 7 28684 -448 89953 
26 2) 32631 435 70972 
ZT 4 317895 -418 63360 
28 3 44210 -482 63928 
28 3 51580 - 460 71180 
30 3 58947 -444 T8517 
aL 3 69474 o42d 88996 
52 3 78947 S45 98289 

Totals: 4396861 5429976 


19.0% savings 
fuel saved in 1000 hrs operation 
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During maneuvering and high speed operations the gas 
turbines would be operated alone with the waste heat 
recovery unit dry. Because of the relatively low gas 
turbine exhaust temperatures involved, it is not expected 


that damage would occur in the WHRU when operated dry. 
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IV. RECOMMENDATIONS FOR FURTHER RESEARCH 


A. WHRU DESIGN MODEL IMPROVEMENT AND EXPANSION 

The type of control imposed on the WHRU design model of 
this thesis was essentially the adjustment of the steam flow 
rate to maintain a superheater steam outlet temperature of 
650 F at the pressure specified. This steam flow rate 
control was further extended to accommodate the maintenance 
Se a Minimum pinch point AT of 25 F. Further investigation 
of WHRU design should include the effects of allowing the 
superheater outlet steam temperature and/or pressure to 
"float" with changing gas turbine exhaust conditions. The 
combination of controls on the superheater outlet pressure, 
temperature, and flow rate could lead to improved WHRU output 
while still maintaining an acceptable minimum pinch point 
temperature difference. 

In order to increase the precision of the model and lead 
to increased knowledge of the gas temperature distribution 
in the WHRU, the present model could be modified for a "fine 
Mesh"approach. That is, instead of solving one pass at a 
time, each pass could be divided into a number of segments. 
The solution could then proceed using either the finite 
difference or finite element method. 

Other possible model improvements are as follows: 

l. Include the calculation of fluid/steam-side pressure 
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2. Additional fin-tube configurations should be considered 
for possible improvement of the heat transfer characteristics 
and the weight and space requirements of the WHRU. 

3. Consideration of inside and outside heat transfer sur- 
face fouling should be given when investigating possible 
fin-tube configurations. 

4. Working fluids other than water should be investigated 
for possible enhanced thermodynamic characteristics at a 


minimum cost in system maintainability. 


B. COGAS MODEL IMPROVEMENTS 

After the waste heat recovery unit, the steam turbine 
and condenser should have the highest priority for modeling. 
Both of these units were modeled as "black boxes" in this 
model. The modeling or performance mapping of "State of the 
art" steam turbines for application in the COGAS model would 
allow the designer to further refine the estimate of system 
performance for various WHRU outputs of steam temperature, 
pressure, and flow rate. The inclusion of a steam turbine 
model would also allow the designer to make size and weight 
estimates for that component. The same benefits would be 
derived from the inclusion of a condenser model in the COGAS 
system. As the condenser pressure is reduced, the system 
performance should improve. The improvement must, however, 
be weighed against the increased condenser size necessary 


at the reduced pressure for constant steam and cooling water 
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conditions. A condenser model would allow the designer 
to make these comparisons. 

In the model of this thesis, pumping power for the 
condensate and feedwater was assumed to be negligible. [In 
the case of the condensate pump, this is probably not a bad 
assumption since it is likely that electric pumps could be 
used, and the cost to the ship's electrical output would 
not be large. The feedwater pump, however, would probably 
be steam turbine driven at some cost to the COGAS system 
output. The required feed pump pumping power for various 
loads should be mapped and included in the model. 

Finally, the likely engineroom machinery layout for the 
COGAS system should be described in enough detail that the 
lengths of piping runs, number of turns, and location of 
valves could be estimated. The inclusion of this feature 
in the model would allow the designer to address pressure 
drop and heat losses between components. This feature would 
also contribute to the COGAS system weight and space predic- 


e10Nn. 


Ce. SYSTEM OPTIMIZATION 

Since the COGAS system is operated with a relatively 
small heat source, optimization should be attempted in the 
design of the system. The technique of non-linear programming 
could be applied to a COGAS system optimization where an 
important system output parameter, say specific fuel con- 


sumption, is expressed as a function of the system design 
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variables such as component weight, component dimensions, 
and gas turbine backpressure. This function is called the 
objective function, and it is maximized or minimized subject 
to constraints which are also expressed as functions of the 
design variables. These constraint functions may be either 
linear or non-linear and express some parameter limit not 

to be exceeded, such as WHRU or condenser total volume. 

For a system of the complexity of the COGAS system a 
"local" optimization approach would probably have to be 
adopted. That is, each of the major sub-models (WHRU, 
steam turbine, and condenser) could be optimized using local 
constraints. Then the entire system could be optimized using 
linking variables, such as enthalpy, flow rates, and tempera- 
tures between components. The process would be repeated 
as Many times as necessary to achieve the final system 
optimization. One available computer-based system using 
the technique of non-linear programming is the COPES/CONMIN 


program [Ref. 14]. 
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APPENDIX B 


This appendix contains the program listings for the 
subroutines called by the WHRU design program and the COGAS 
model to obtain steam, water, and air properties. Refer- 
ences 15 and 16 provide relationships for use in calcula- 
ting the thermodynamic properties of steam and water. The 
relationships of Ref. 15 for the thermodynamic properties 
of steam were used in Subroutine SS of this appendix. 

The remainder of the subroutines for steam, water and 
air properties were obtained by fitting polynomials to the 
values for those properties found in References 17 and 18. 
In all cases, the polynomial coefficients were obtained 
using a "canned" program for least squares regression pro- 
vided with the Hewlett-Packard System 9845 Computer installed 
in the Naval Postgraduate School Mechanical Engineering 


Department. 
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